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The tube-like cages of stiff polymers in entangled solutions have been shown to exhibit charac- 
teristic spatial heterogeneities. We explain these observations by a systematic theory generalizing 
previous work by D. Morse (Phys. Rev. E 63:031502, 2001). With a local version of the binary 
collision approximation (BC A) , the distribution of confinement strengths is calculated, and the mag- 
nitude and the distribution function of tube radius fluctuations are predicted. Our main result is a 
unique scaling function for the tube radius distribution, in good agreement with experimental and 
simulation data. 
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I. INTRODUCTION 

Entangled solutions of stiff polymers are minimal 
model systems to generate a fundamental understand- 
ing of the origin of the mechanical properties of the cy- 
toskeleton. This complex polymer scaffold maintains the 
stability and integrity of animal cells and is comprised of 
three types of semiflexible filaments, microtubules, actin, 
and intermediate filaments, with backbone diameters in 
the nanometer range but persistence lengths on the order 
of 10 _1 . . . 10 3 fim [TH5]. Single stiff biopolymers exhibit a 
rich mechanical response [JHZ] • In- vitro reconstituted so- 
lutions of such biopolymers hint at how cells can acquire a 
considerable macroscopic strength from a purely topolog- 
ical microscopic constraint and thermal fluctuations, uti- 
lizing a minimum amount of material. Though the indi- 
vidual polymers only have to respect a simple constraint, 
namely the mutual impenetrability of the polymer back- 
bones, complex soft-solid mechanical behavior arises at 
densities that would correspond to a very dilute gas with- 
out polymerization and a certain flexibility allowing for 
thermal backbone undulations. To deform an entangled 
polymer, surrounding polymers need to be pushed out 
of the way, as familiar from knotted strings. This mech- 
anism leads to confinement of the individual polymers 
in effective tube-like cages, from which they only escape 
very slowly by a snake-like motion called reptation [5] • 
The suppression of chain motion perpendicular to the 
tube backbone is responsible for the remarkable integrity 
of the transient network. A microscopic derivation of 
this confinement poses formidable theoretical challenges, 
and there has so far been little progress beyond the in- 
troduction of basic topological invariants characterizing 
polymer entanglement [101 111] and a phenomenological 
primitive path analysis |12j . 

Nevertheless, self-consistent approximations for the 
dynamics of rigid |T3l [14] and the equilibrium statisti- 
cal mechanics of semiflexible [TS] topologically entangled 
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FIG. 1. (Color online) Test polymer in a background solution, 
confined into a tube of spatially varying radius R(s). The 
chemical distance L e indicates the characteristic scale of the 
tube heterogeneities. 



chains have been worked out, and these treatments can 
predict salient properties of the reptation dynamics and 
the postulated tube. For stiff but not rigid polymers, 
with a large but finite persistence length l p , the tube 
confinement of the transverse fluctuations of a represen- 
tative test chain is implemented by a harmonic potential 
of stiffness <fi. The confinement geometry is character- 
ized by an entanglement (or collision) length L e <C l p 
and the tube width R L e [TH1 [T7] , both of which are 
functions of </>. The former is a measure of the aver- 
age spacing between adjacent collisions with background 
polymers, each of which contribute an amount fegT to 
the average confinement energy of the test chain. The 
latter measures the magnitude of the confined thermal 
fluctuations (Fig. [I]) . Its mean value R as a function of 
monomer concentration has been predicted theoretically 
based on a binary collision approximation (BCA) and an 
effective medium approximation (EMA) [15] . The BCA 
focuses on the pairwise entanglement topology of a test 
chain, while the EMA aims to account for the collective 
network fluctuations. More recently, these mean-field 
type theories have been challenged by the observation of 
pronounced heterogeneities of the local tube width R(s) 
along the tube contour, which have been systematically 
studied in experiments [THHH] and in simulations |23j . 
The tube heterogeneities have been statistically quanti- 
fied by a broad and skewed tube width distribution P(R), 
which has been analyzed by an empirical model |21j and 
by a generalization of the BCA [22] . 
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In the following, we develop a systematic, BCA-based 
theory to describe the fluctuations of the tube radius 
in an entangled polymer solution on the scale of indi- 
vidual tube collisions. Thereby, the local tube radius 
heterogeneities R(s) and their distribution P(R) can be 
determined, and P(R) is found to be a universal non- 
Gaussian scaling function with a stretched tail. By com- 
parison with the segment fluid approach of Ref. |22j . in 
which the entangled solution is effectively mapped onto 
an ensemble of entanglement segments, we predict the 
segment length L (which was previously treated as a fit 
parameter). The magnitude of the tube width fluctua- 
tions is compared with published experimental data. In 
Ref. [21] . P(R) was instead estimated based on an ad-hoc 
distribution of the local mesh size. The result turned out 
to be unphysical at small values of i?, however. As we 
show in the following, the fluctuations of the tube radius 
R can be comprehensively described without additional 
assumptions based on a generalization of the BCA. 

The remainder of the paper is organized as follows. 
In Sec. [Tl] we introduce the fundamental concept of the 
tube and the wormlike chain (WLC) model for a single 
confined scmiflexible polymer. In Sec. |III[ we then sum- 
marize the basic assumptions underlying the BCA as an 
approximation to the topological problem. Subsequently, 
in Sec. |IV[ we discuss the statistical distribution of the 
confinement strength, which explains the fluctuations of 
the tube radius, that are derived in Sec. [V] In Sec. |VI[ 
the magnitude of tube radius fluctuations is used as an 
input to the segment fluid model, which predicts a scal- 
ing function for the tube radius distribution P(R). The 
analytical results are compared to experimental data in 
Sec. ECO 



II. BASIC ELEMENTS OF THE TUBE MODEL 

A. Time scale separation and topology 

We consider a stiff test polymer in the presence of sur- 
rounding uncrossable polymers, which are imposing topo- 
logical constraints on its conformation. We restrict our 
discussion to tightly entangled polymers that are charac- 
terized by small transverse excursions around an average 
path, the preferred contour. The tube concept concerns 
quantities in an intermediate equilibrium, i.e. on time 
scales t £ « i « tj, where r e is the time for the confined 
degrees of freedom to equilibrate inside the tube and Td 
is the disengagement time of the polymer from its initial 
tube. In what follows, a strong scale separation r e <C Td is 
assumed. Then, in the idealized limit Td — > oo, the topo- 
logical relationships of the solution will be asymptotically 
conserved, and the average positions of the background 
polymers and their mutual topological relationships can 
be considered as effectively frozen (or "quenched" ) , thus 
collectively giving rise to a (quasi-static) confinement po- 
tential representing the tube. We denote the thermal av- 
erage with respect to a given "quenched" configuration 



by angular brackets (■ ■ ■ ) and the average over different 
configurations and topologies of the tube by an over-bar 
~. These ensemble averages correspond to temporal av- 
erages over several time intervals of length r e and Td, 
respectively. 



B. Statistical mechanics of a single entangled stiff 
polymer 

1. Transverse distance distribution 

To describe the physical properties of the test polymer, 
we assume that the effect of confinement can, to leading 
order, be described by a harmonic confinement potential. 
Hence, we use the weakly-bending Hamiltonian 
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for the transverse fluctuations rj_(s) of the test polymer 
about the straight ground state of a rigid rod, with a 
local confinement strength (f>(s) that will be determined 
self-consistently. We use natural energy units (/csT = I), 
such that the persistence length l p = n/ksT is synony- 
mous with the bending rigidity n. We define the arc- 
length dependent tube radius R(s) via the variance of 
one component of the confined transverse fluctuations, 
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Approximating the free energy by an effective Hamil- 
tonian Hconf that is quadratic in the fluctuations is 
equivalent to approximating the distribution of P[rj_(s)] 
in a given configuration by a Gaussian. Experiments 
QH[2T][22] and simulations [HIM] indicate that the dis- 
tribution of transverse distances is indeed Gaussian for 
small transverse displacements rj_(s) on the order of the 
tube radius R. It can be shown theoretically, that this 
assumption is in accord with a self-consistent treatment 
of the tube [15] . 



2. Tube radius R and entanglement length L e 

As a first step, we consider the case of a test polymer 
in a homogeneous (cylindrical) tube that can be charac- 
terized by the spatial average of a local confinement 
strength </>(s), and hence a tube radius R(s) = Rq = 
const. The variance of rj_(s) [Eq. |2])] is obtained from 
the tube Hamiltonian Eq. ([I]) via equipartition, such that 

R 2 (s)=[P^ (3) 
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is the square of the tube radius corresponding to a homo- 
geneous confinement strength <j>. Heterogeneities of the 
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FIG. 2. (Color online) Topology of two confined polymers 
inside their tubes at a collision point: for fixed preferred tube 
contours in the transparent state, the polymers are found ei- 
ther in the 'above' (+) or 'below' (— ) configuration. Adapted 
from Ref. [15]. 

tube potential and of the tube radius are discussed below 
(in Sections |IV| & [V} , where we show how small spatial 
fluctuations 6<f>(s) = <f)(s) — 4> lead to spatial variations 
of R(s) about its average value R. We assume in what 
follows that the peak of the corresponding distribution 
P(R) is sufflcently well defined such that the average R 
and the typical value Rq can be used interchangeably. 

The second characteristic quantity of the tube geome- 
try, the entanglement length L e , is defined by assigning 
a harmonic confinement energy equal to kgT (= 1 in 
our units) to every collision, and identifying L e with the 
collision length. Writing <f> for the average confinement 
potential strength in Eq. ([!]), equipartition yields 




III. BCA 

We recapitulate the essential arguments that are 
needed to derive the BCA and to understand the rea- 
soning that follows. 

The BCA was designed as an approximation to the un- 
derlying topological many-body problem, suitable for es- 
timating the absolute value of the average tube radius R, 
self-consistently. One considers an elementary encounter 
('binary collision') between two tubes, calculates the free 
energy of confinement due to the uncrossability of the 
chains, and sums over possible configurations of the pair 
of tubes. The key approximation consists in neglecting 
correlations between multiple collisions. 

The dynamic entanglement problem is cast into equi- 
librium statistical mechanics language by assuming that 
the tube contours are (temporarily) frozen. In such a con- 
figuration, an arbitrarily chosen point on the tube con- 
tour associated with the test chain is characterized by a 
Gaussian distribution of transverse distances, as outlined 
in the preceding section. Its standard deviation is ap- 
proximated by the average tube radius R. Consider now 



a second background chain passing within a distance R 
from the chosen point on the test chain. We refer to this 
event as a tube collision. To calculate the contribution 
of the pair collision to the confinement free energy, the 
BCA distinguishes between two states: a hypothetical 
state, in which the test chain is transparent with respect 
to collisions with the background chain and a state in 
which the chains are mutually uncrossable. Due to the 
key assumption that the collisions along the test chain 
are uncorrelated, the environment of the collision point is 
completely random in the transparent state, i.e. the dis- 
tribution of transverse distances is unchanged from the 
Gaussian distribution with standard deviation R. The 
average free energy of confinement follows from a topo- 
logical argument. In the uncrossable state for this pair of 
chains, the configuration space of transverse fluctuations 
is divided into two disconnected subspaces correspond- 
ing to the 'above' (+) and the 'below' (— ) configuration, 
as depicted in Fig. [2] The average free energy in the 
uncrossable state is therefore obtained by averaging the 
confinement potential in each of these subspaces over the 
probability for a specific topology and collision geometry. 
The resulting average free energy is equated with that of 
the transparent state to obtain a self-consistent estimate 
of the tube radius R. 

For the original mathematical implementation of the 
above ideas we refer the reader to Ref. [TS]. We note 
that a slightly corrected estimate for the prefactor in the 
scaling result for the average tube radius R with concen- 
tration was calculated in Ref. [53] (supplement). 



IV. DISTRIBUTION P[<f>] OF TUBE 
STRENGTHS 

The original BCA is exclusively concerned with aver- 
age values R, (p. The scaling of these values with con- 
centration p has been tested experimentally [501 H3 12S] , 
but the prefactor is sensitive to the precise control of the 
experimental conditions and is usually treated as a fit 
factor. For a more detailed comparison of theory, experi- 
ment and simulations, knowledge not only of the average 
value but of the richer and more robust tube radius dis- 
tribution P(R) is desirable (cf. Sec. |V|. 

As a first step towards calculating P(R), we derive the 
distribution P[(f>] of the local confinement strength <f>(s) 
of a test polymer. Morse [TS] gave an explicit expression 
for the confinement free energy of a test chain collid- 
ing with a medium chain. We will explicitly adopt the 
mathematical approximation of straight tube contours, 
that was implicitly made in the previous work, and it 
is shown that inconsistencies resulting from this approxi- 
mation are avoided by verifying that the calculated quan- 
tities do not depend on the overall chain length. Let the 
distance of shortest approach between two preferred con- 
tours with orientations u, u' and centers of mass r, r' be 
x = (r — r') -e x , where e x = ux u'/|ux u'| is the direction 
perpendicular to both preferred contours, then the free 
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FIG. 3. (Color online) Overlap area (parallelogram) of two 
tubes of length L, represented by their preferred contours with 
orientations u, u' enclosing an arbitrary angle, and charac- 
teristic function \. 



energy of the test polymer whose preferred contour has 
been uniformly displaced by a vector h = he^ is 



F± (h) = - In $ ± 



x — h cos tp 
R 



(6) 



Here, the sign ± refers to the specific topology (cf. Fig. [2]) 
and cosV' = e/, • e x . The function §(x) is given by the 
restricted partition sum of the Gaussian fluctuations of 
the test polymer in the presence of an uncrossable test 
chain ITBI. 
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It can be interpreted as the probability p±(x) = 
&(±x/R) of finding a specific topology. 

From Eq. @, we obtain the confinement strength <f> in 
a given configuration of preferred contours and topology 
as the second derivative of the free energy, 
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To derive the distribution of confinement strengths 
4>{s), we now turn back to the central BCA approxima- 
tion that the collisions between the test polymer and the 
background polymers are independent localized events. 
Due to the requirement R <C L ei we may, without loss of 
generality, even treat them as point-like and express the 
confinement potential <fi(s) per unit length at a point s 
on the test polymer as 



N 

<t>(s) = ^x(r i ,u l ,u)(5(s - z l )4>±{x l ). 

i=l 



(9) 



Here, %(rj, Uj, u) is a characteristic function of overlap 
between the colliding tubes, which takes on the value 
one whenever a tube collision occurs, and zero otherwise 



(for a graphical definition see Fig. (31). We introduce it, 
here, merely as a convenient tool to facilitate the formal 
manipulation of the following expressions. The coordi- 
nate Zi in the argument of the (5-function is the point of 
shortest approach (the collision point) between the two 
tubes on the test polymer. 

The distribution P[<f>] of confinement strengths that 
follows from Eq. ^ is of the Holtsmark type [2H |57] and 
describes the total confinement potential of the test chain 
as a sum of contributions resulting from uncorrelated 
collisions with medium chains. Explicitly, its moment- 
generating functional is given by (cf. App. [A]) 



P[w(s)] = exp 



X e tw[z )0±(s) p±(a/) 
(10) 

From this, we obtain the average of 4>(s) by functional 
differentiation of the logarithm of Eq. ( 10 ) with respect 
to the field w(s) (cf. App.lA), 
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The integral in Eq. (Ill is numerically evaluated and 
gives 4> = anL/R, with a — 0.502, in agreement with 
earlier results |15j (applying the slight numerical correc- 
tion discussed in Ref. [22 , supplement). 

We proceed analogously to obtain the second cumulant 
(cf. App. [A|, the correlation function 



(s)<l>(s>)=p~5( S -s'), 
R 
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where (3 — 0.0941. The uncorrelated character of the 
collisions is apparent from the 5-function on the RHS 



of Eq. (12). Both should be understood in the coarse- 
grained sense, assuming <j)(s)(f>(s') = when s — s' 3> R. 



V. GAUSSIAN APPROXIMATION TO THE 
TUBE RADIUS DISTRIBUTION P{R) 

We can now turn the distribution P[4>], which we char- 
acterized by its first two cumulants, into a Gaussian ap- 
proximation to the tube radius distribution P(R), by 
calculating the linear response of the local tube radius 
R(s) to spatial changes (heterogeneities) in <p(s). We 
begin with the observation that the correlation func- 
tion of the (projected) transverse fluctuations, C(s, s') — 
(rj_(s)-rj_(s'))/2 obeys the following differential equation 
(cf. App. Ill 



l p d^C(s, s') - <^(a')C(a, s') = S(s - s') 



(13) 



The tube radius is given by R(s) = [C(s, s)] 1 / 2 and the 
linear response expression for this quantity is calculated 
in App. [B] as 



R(s) =R 



1 

2! 



ds'G 2 (s-s')5<f>(s'). (14) 



5 



The variance of the tube radius R(s) at a randomly 
chosen point s on the test polymer is now calculated from 



Eq. (14) as 
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Within the BCA, with its trivial spatial correlations of 
the confinement potential </>(s), Eq. (12), this reduces to 



5R 2 = p 
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The integral Eq. ( |16[ ) is numerically evaluated using an 
explicit expression for G(s) [Eq. (B5|], 



dsG 4 (s) = 7 i? : L e , 



(17) 



where 7 = 0.5125, which yields the final result for the 
variance of P(R), 
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Using the self-consistent solutions for R = 

(4a)" 3 / 5 (nL)- 3 / 5 Zp 1/5 and L e = (a/8)~ 2 / 5 (nL)~ 2 / 5 ll /5 
|15j . mean and variance of the Gaussian approximation 
to the tube radius distribution P(R) arc completely 
determined in terms of the contour length concentration 
nL and the persistence length l p . In particular, the 
coefficient of variation c v = [(5R) 2 ] 1 / 2 / R turns out to be 
a concentration-independent constant, 




= 0.155. 



(19) 



VI. SEGMENT FLUID APPROXIMATION 

In Ref. [22], a broad distribution P(R) of the tube 
radius was found. The derivation of an analytical re- 
sult for this distribution was also based on a Holtsmark- 
type distribution for the confinement strength <f> resulting 
from uncorrelated collisions, but the latter were averaged 
over the characteristic length L of entanglement segments 
(which is why the approach was called a "segment-fluid" 
approximation). It was argued that this length is on the 
order of L e . We now show that this choice is indeed jus- 
tified and predict the precise value of the segment length. 

The tube radius distribution P(R) was given as an 
analytical approximation in Ref. [22], 
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where k = 4.013n' L 2 R, n' is the number density of en- 
tanglement segments and T(x) is the Gamma function. 




R/R 

FIG. 4. (Color online) Reduced distribution RP(R/R) of the 
tube radius R. The solid line is the segment-fluid prediction, 
the dashed line is a Gaussian approximation with coefficient of 
variation c v given by Eq. (19 1. Inset: Coefficient of variation 



c' v of the segment fluid approximation versus reduced segment 
length e = L/L e , and the predicted value. 



It was noted [22] that P(R) can be written as a scal- 
ing function P(R) = (l/R)f(R/R,n'L 2 R). This implies 
that the coefficient of variation is given by a constant, 
c', = g(n'L 2 R). If we make the ansatz L — eL e , with 
an undetermined dimensionless constant e and use the 
self-consistent values for R and L e given above, we get 
c' v = g(e/a). The function g(x) is easily evaluated nu- 
merically and is shown in Fig. [4] (inset). To fix e and 
thus the length L of entanglement segments, we require 
cL = c v and obtain numerically e = 0.85. 



Evaluating the tube radius distribution P{R) [Eq. ( 20 )] 
at this value of the reduced segment length L/L e , we 
obtain k — y(R/R) s / 3 — 6.79 and hence all parameters 
occurring in P(R) are now fully specified, giving 



P{R) = Nexp 
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(21) 



with a normalization constant N — 2.434 x 10 3 . The cor- 
responding unique reduced distribution RP(R) is shown 
in Fig.[4]and compared to the Gaussian distribution with 
the same c v . Beyond what was achieved in Ref. [22], the 
functional form of P(R) is now fully determined. As 
can be seen in Fig. |4j the distribution P(R) is positively 
skewed and has a broad tail at large values of R. 



VII. COMPARISON TO EXPERIMENT 

The functional form of P(R) was compared to experi- 
mental data in Ref. [22J , and very good qualitative agree- 
ment was found, using the value of the segment length L 
as a fit parameter. Remarkably, also our above predic- 
tion of a constant value for the coefficient of variation c v , 
which can be checked against a whole set of independent 
measurements, is nicely confirmed by the data (Fig. [5]). 
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tion, for which we gave an analytical approximation in 
Eq. (21 1. Our results provide a quantitative characteriza- 
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FIG. 5. (Color online) Comparison of tube radius fluctua- 
tions in different experiments and simulations of semidilute 
solution s of F-actin. Shown is the coefficient of variation 
c v = V 8R 2 /R for the fluctuations of R. Symbols corre- 
spond to data taken from the literature. Circles: experimental 
data taken from Ref. 22 (Fig. 3, inset); squares: experimen- 
tal data from Ref. [28]; diamonds: experimental data from 
Ref. pjj] for the fluctuations of the response coefficient a±, 
converted to fluctuations of the tube radius R; upright trian- 
gles: experimental data from |21j (Fig. 2); downward facing 
triangles: simulation data from Ref. |23| (Fig. 9). The solid 
dashed line is the prediction of Eq. 
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The figure summarizes literature data for c v against 
monomer concentration c from various experiments and 
simulations of semidilute actin solutions. The dashed 
line is our prediction from Eq. (19). Two results are 



evident from this plot. First, the data scatter within a 
band of c v — 0.2-0.4. Second, the theoretical prediction 
lies below most of the data points and thus provides a 
lower bound for the observed tube radius fluctuations. 
This suggests that a constant value for the coefficient 
of variation is indeed consistent with the reported data, 
but that the heterogeneities are actually about twice as 
strong as predicted. 

This is not entirely unexpected, since the BCA, on 
which our theoretical derivation relies, is not meant to 
describe the absolute value of the tube radius quanti- 
tatively. In fact, the BCA alone is well known to un- 
derestimate the tube fluctuations, since it does not take 
into account the collective fluctuations of the surround- 
ing medium into which the tube is embedded |15j . Cor- 
responding quantitative discrepancies with experiments 
have been reported before [2"5] . 



VIII. CONCLUSION 

We have calculated the fluctuations of the tube radius 
in entangled solutions of semiflexible polymers, based on 
the binary collision approximation (BCA). We predict 
that the shape of the tube radius distribution is given 
by a universal (concentration-independent) scaling func- 



tion of the local packing structure of entangled biopoly- 
mer solutions in terms of distribution functions, which 
are at the same time a more sensitive and more robust 
means for comparing data and theory than average values 
alone. We hope that the methods of analysis established 
here may find application in future experimental studies, 
e.g. in microrheology 29, 30 , or in the interpretation of 
simulation data [23j [31] . Further theoretical questions, 
such as the characterization of the distribution of tube 
contours [201 122] > are currently under study. 
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Appendix A: Calculation of the distribution of 
confinement strengths P[<f>] 

The distribution of the local confinement strength 
formally obtained as the average P[<p(s)] = 



5[4>(s) — <p(s)] over all possible configurations of tube con- 
tours and topologies. The corresponding characteristic 
functional, P[w(s)], follows from a functional Fourier 
transform and Eq. Q as 



P[w(s)] = exp 




(Al) 



Here, tangents to the test and the background chains' 
tube backbones are denoted by u and u i; and the vec- 
tor Ti connects the tubes' center of mass of the test 
chain with that of the background chain i. Its coor- 
dinates Xi, yi and Zi are defined along the directions 
e x .i = (ux Uj)/|ux Uj|, e Vii = uxe X)i and e Xii = u. The 
quenched average 7TT is implemented as the simultaneous 
average over the probability p±(xi) = <&(±Xj) of finding 
a specific topology '+' or '-' [defined after Eq. Q] and 
over the uniformly distributed centers of mass and orien- 
tations of the tubes. Since the chains' preferred contours 
are assumed to be uncorrelated, the average over the N 
preferred contours and the topology of the background 
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chains relative to the test chain factorizes as 
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Here, the integral over orientations extends over the half- 
sphere (indicated by a prime) . Exploiting the formal def- 
inition of x as a characteristic function of overlap, which 
amounts to setting the factor in the brackets for non- 
overlapping chains to unity (since the probability p±(x') 
is normalized), we get 



P[w(s)} 



E 

± 



dr' 
~V 



' da' 
27 



iV 



(A3) 

Using n — N/V for the polymer number concentration 
and performing the limit N — > oo, Eq. (10) in the main 
text is obtained. 

The first cumulant is obtained by functional differen- 
tiation of the characteristic functional Eq. ( 10 ) with re- 
spect to the field w(s), 



8w(s) 
nL % 
2n ^ 



In P[w{s)} 



(A4) 



w(s)=0 

du' I dr' 5(z' - s) X p±(x')<j> ± {x'). 

(A5) 



Using the fact that the integral of x${ z ' ~ s ) over dy' 
and dz' is the height LsinO of the overlap area (Fig. [3]), 
carrying out the second derivative of — In $ and the an- 
gular integrals, and using ^ ± $(x') = 1, one arrives at 
Eq. (11) in the main text. Analogously, we obtain the 



second cumulant, the correlation function 



5w(s)Sw(s') 
nL 
2^ 



-lni>(a)] 



w(s)=0 



(A6) 
(A7) 



^ J du' J dr' 8(z' - s)5{z' - s') X P±(x')tf ± {x'). 

(A8) 



Applying a similar reasoning as above to simplify the 
equation and numerically evaluating the remaining x- 
integral, we obtain Eq. (12) in the main text. 



Appendix B: Heterogeneous tube radius R(s) 



— In Z[f±(s)} of a confined WLC in the presence of an ex- 
ternal transverse force fj_(s). The corresponding Hamil- 
tonian is H = H con{ + J ds f±(s)r±(s). Since 



(r ± (s)r ± (s')) 



8\\\Z 

Wis)' 
8\nZ 



8f ± (s)8f x {s') 
S(r ± (s)) 



f±{s)=0 



fx(s)=0 



(Bl) 
(B2) 
(B3) 



it follows that C(s,s') = (r±(s)r±(s')) is the functional 
inverse of C(s, s') _1 = 8i±(s)/8(r±(s')) . Since the force 
fj_(s) = — (5H con f/5r±(s)) producing an average dis- 
placement (rj_(s)) is given by —l p dg(r±(s))~(j)(s)(r±(s)), 
Eq. ( 13 ) follows by partial integration. 



A solution of Eq. ( 13 ) would exactly describe the tube 



heterogeneities that follow from a heterogeneous con- 
finement potential (f>(s). However, no such solution is 
available for arbitrary <j)(s). Therefore, we write 4>(s) = 
(j) + 6<f>(s) with small fluctuations 8<fi(s) about the aver- 
age confinement strength <f>. A simple first-order pertur- 
bation scheme for C(s, s') = C^(s, s') + SC{s 7 s') is set 
up by requiring C*- -* to be the response function in the 
homogeneous case, where 5<p(s) = 0, 



C^(s-s')=G(s-s')= J 



dq e tq(s - s '^ 
2tt Lg 4 + 0' 



(B4) 



The explicit expression for the Fourier transform in 
Eq. (B4) can be obtained analytically and written [us- 
ing Eqs. @ & <(5j|] as 



G(S-S') =i?2 e -2| S |/ io 



cos 2 



sin 2 



(B5) 



The leading-order response 6C(s, s') to the perturba- 
tion 6<f>(s, s') is obtained from Eq. ( 13 ) if small terms 
0(6C8<j>) are neglected, 



5C(s,s') 



ds"G{s - s")S<l>(s")G(s' - s"). (B6) 



The fluctuation-response relation Eq. ( 13 ) for the cor- 



relation function C(s, s') is derived from the free energy 



Eq. ( 14 ) in the main text is obtained by writing R(s) 
R + 8R(s) with 5R(s) = S(R 2 )/2R = 8C(s lS )/2R . 
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